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^ \ Abstract 

O' 

We introduce a semi-parametric estimator of the Poisson intensity parame- 
ter of a spatial stationary Gibbs point process. Under very mild assumptions 
satisfied by a large class of Gibbs models, we establish its strong consistency 
1 and asymptotic normality. We also consider its finite-sample properties in a 

X/} • simulation study. 
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1 Introduction 

Spatial Gibbs point processes form a major class of stochastic processes allowing the 
Tj- 1 modelling of dependence of spatial point patterns. They have applications in many 

O ■ scientific fields such as biology, epidemiology, geography, astrophysics, physics, and 

economics. General references covering as well theoretical as practical aspects of 
these processes are e.g. [20], [H] or [T2] . 

Gibbs point processes in M. d can be defined and characterized through the Pa- 
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pangelou conditional intensity (see [H]) which is an application A : R x Q — y 
where Q is the space of locally finite configurations of points in M. d . The Papan- 
gelou conditional intensity can be interpreted as follows: for any u G M. d and x e Q, 
X(u, x)du corresponds to the conditional probability of observing a point in a ball 
of volume du around u given the rest of the point process is x. In the present paper, 
we assume that the Papangelou conditional intensity is decomposed as 

A(u,x) = /3A(u,x), (1.1) 

where /3 is a positive real parameter and where A : M. d x Q — y R + . If the point 
process corresponds to a homogeneous Poisson point process A(tt, x) = 1. Therefore, 
we suggest the name Poisson intensity parameter for the parameter 0. The goal 
of this paper is to develop an estimator of without specifying anything on the 
function A. Our semi-parametric estimator is based on a smart application of the 



1 



Georgii-Nguyen-Zessin formula (see j9] for a general presentation and Section 12. ip 
which is also a way to characterize a Gibbs point process. We propose a very simple 
ratio-estimator of /3 based on a single observation of a stationary spatial Gibbs 
point process observed in a bounded window and we prove its asymptotic properties 
(strong consistency and asymptotic normality) as the window expands to M. d . The 
problem considered in this paper is related to two different works that we detail 
below. 

First, due to its easy interpretability the most popular Gibbs model is without 
any doubt the (isotropic) pairwise interaction point process. Its Papangelou condi- 
tional intensity is given by 

X(u, x) = p\(u, x) with \(u,x) = Y[g(\\v -u\\). (1.2) 

t)€x 

The function g is called the (isotropic) pairwise interaction function. Estimating the 
model (ll.2p consists in estimating (3 and the function g. To the best of our knowledge, 
this awkward problem has been considered only in [8]. The authors propose several 
methods to estimate non parametrically the function g. Our objective is to make 
some advances in that direction which justifies first to propose an estimator of /3 
independently of the function g and to understand its properties. 

Second, many popular methods exist to estimate a parametric Gibbs model. This 
includes methods based on the likelihood or the pseudo-likelihood and the Takacs- 
Fiksel method. We refer for instance to [14] . [T2] or [13] and the references therein for 
an overview of these methods. A recent contribution [1] investigates a new approach. 
The following exponential model is considered 

\(u, x) = /3\(u, x) with log \(u, x) = 9 T v(u, x) (1.3) 

where 9 is a real p-dimensional parameter vector, v(u, x) = (v i(u, x), . . . , v p (u, x)) T 
for measurable functions — > R, i — 1, . . . ,p. PQ develops a new estimator 
of the parameter 9 where the parameter /3 is treated as a nuisance parameter. The 
main advantage of the estimator is that it is obtained via the resolution of a linear 
system of equations whereas all previous methods cited above require an optimiza- 
tion procedure. The ratio-estimator of we propose is also very computationally 
cheap. Hence, it could serve as a complement of the estimate of 9 proposed in [TJ to 
completely identify the model (11.30 . 

The rest of the paper is organized as follows. Section [2] gives the background and 
notation, presents examples of Gibbs point processes and our main assumptions. 
Section [3] deals with the core of the paper. We present our estimator and derive its 
asymptotic properties as the window of observation expands to M. d . A simulation 
study is conducted in Section H] where we illustrate the efficiency of our estimator. 
Proofs of the results are postponed to Appendix. 
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2 Background and assumptions 



2.1 Basic definitions and tools 

A point process X in M. d is a locally finite random subset of IR d , i.e. the number 
of points N(A) = n(K\) of the restriction of X to A is a finite random variable 
whenever A C M d is a bounded region (see We let Q be the space of locally 
finite point configurations. If the distribution of X is translation invariant, we say 
that X is stationary. 

The Papangelou conditional intensity completely characterizes the Gibbs point 
process in terms of the Georgii-Nguyen-Zessin (GNZ) Formula (see [T7] and for 
historical comments and [9J or [TJ] for a general presentation). The GNZ formula 
states that for any measurable function h : M. d x fl such that the left or right hand 
side exists 

E V/i(u,X\u) = E / h{u,X)\{u,X)du. (2.1) 

We will not discuss how to consistently specify the Papangelou conditional intensity 
to ensure the existence of a Gibbs point process on but rather we simply assume 
we are given a well-defined Gibbs point process. The reader interested in a deeper 
presentation of Gibbs point processes and the existence problem is referred to [T9j [T8] 
or [7], see also Section l2~2l for a few examples. 

The concept of innovation for spatial point processes is proposed in [2] and is a 
key- indredient in this paper. Inspired by the GNZ formula, it is defined as follows: 
for some measurable function h : M. d x Q — > M, the h— weighted innovation computed 
in a bounded domain A is the centered random variable defined by 

J A (X,/i)= h{u,X\u)- I h{u,X)X{u,X)du. (2.2) 
«ex A 

We end this paragraph with a few notation. The volume of a bounded Borel set A 
of R d is denoted by | A |. The norm II • II stands for the standard Euclidean norm. For 
two subsets Ai and A 2 of M. d , d(Ai, A 2 ) corresponds to the minimal distance between 
these two subsets. Finally, B(u,R) denotes the closed ball centered at u G W 1 with 
radius R. 

2.2 Examples of Gibbs point processes 

We present in this paragraph several classical examples through their Papangelou 
conditional intensity. Further details on these models can be found in [2] for exam- 
ples (i)-(iv) and (vi)-(vii) and to [TT] for the example (v). Let u G x G Q and 
R>0. 

(i) Strauss point process. 

A(w,x) = /3 7 ™[o,*](«> x ) 

where (3 > 0,7 G [0, 1] and ri[ 0iR ](u, x) = yjl(||t; — u|| < R) represents the number 
of i?-closed neighbours of u in x. 
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(ii) Strauss point process with hard-core. Let < 5 < R the Papangelou condi- 
tional intensity is 

f if d(u,x) < S 

A(«,x) - | ^ n]SM{u>x) otherwige 

where /3 > 0, 7 G [0,1]. 

(iii) Piecewise Strauss point process. 

v 

A(u,x) = p J_J_7y 

i=i 

where /3 > 0,7,- G [0, 1], n[ Rj _ liR:j] (u, x) = - u|| G [#j_i, i^-]) and i? = < 

i?X < . . . < i? p = i? < 00. 

(iv) Triplets point process. 

A(w,x) = /3 7 ( s [o,fl]( xUu )- s [o,fl]( x )) 

where (3 > 0,7 G [0,1] and S[o,i?](x) is the number of unordered triplets that are 
closer than R. 

(v) Geyer saturation point process with saturation threshold s 

A(w,x) = /3 7 *( x u«M(x) 

where (3 > 0,7 G [0,1], s G N* and i(x) = ^^ ex min(s, 7i[ 0j fl/2](u, x \ «)). 

(vi) Lennard- Jones model with finite range. This corresponds to the pairwise 
interaction point process defined by 

\(u,x)=Pl[g{\\v-u\\)) 

DgX 

with for r > 0, g(r) = exp (6 6 r- 6 - # 12 r- 12 )l ]0 ,/?](r), 6 > and /3 > 0. 

(vii) Area-interaction point process. 

A(u, x) = /3 7 ^(-u«)-A(x) 

for /3 > 0,7 > and where A(x) = | U„ ex R/2)\. 

Examples (i)-(iii) and (v)-(vi) are pairwise interaction point processes whereas 
example (iv) (resp. (vii)) is based on interactions on cliques of order 3 (resp. of 
any order). Example (vi) is the only example which is not locally stable but only 
Ruelle superstable, see [19]. We recall that the local stability property states that (for 
stationary Gibbs models) the Papangelou conditional intensity is uniformly bounded 
by a constant (see [14]). 

2.3 Main assumptions and discussion 

In this section, we consider the model (11. ip and we present the assumptions we 
require to propose in Section 13.11 an estimator of the Poisson intensity parameter. 
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To stress on the fact that the Papangelou conditional intensity satisfies (11. ip and 
that we aim at estimating f3 we denote, from now on, the corresponding Papangelou 
conditional intensity A^. 

The two following assumptions will be considered and discussed throughout the 
paper. 

(a) The Papangelou conditional intensity is invariant by translation has a finite 
range R, i.e. 

Xp(u,x) = \p(0,T u x) and \p(u, x) = Xp(u, *b(u,r)), (2.3) 

for any u G M d , x G ft, (3 > 0. 

(b) The function A satisfies for any u G M d 

A(u,0) = l. (2.4) 

We emphasize that no more assumption will be required. In particular we do not 
require that the Gibbs point process is locally stable. This is noticeable since local 
stability is often an assumption made when dealing with such models, see e.g. [3J. 
We also underline that in the following, we do not assume that R is known. In place 
of this, we assume we are given an upper-bound of R. 

Let us discuss these assumptions regarding the examples presented in Section l2~2l 
Examples (i)-(vii) exhibit Papangelou conditional intensities which satisfy (II. ip and 
(12. 3p with finite range R. Example (i)-(vi) satisfy the assumption (12.41) . whereas 
example (vii) is a counter-example since for the latter one 
for 7 > 0. 

3 Semi-parametric estimator of f$* 
3.1 Definition 

In the following, we let /3* > denote the true parameter and we denote by Pp* the 
probability measure of X, i.e. X ~ Pp*. 

Our estimate is based on an application of the GNZ formula ( 12. ip for the following 
function h given for all u G R d , x G ft and R > R by 

%,x) = t(d(u,x) >R) = t(*n B(u,R) = 0). (3.1) 

The idea is formulated by the following result. Let A be a bounded region of R d and 
let ]Va(X; R) and Va(X; R) denote the following random variables 

iV A (X; R) = 1((X \ u) D B(u, R) = 0) 

ttgX A 

V A (X;R)= [ l(Xn%H) = 0)dtt. 

J A 
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Proposition 3.1. Assume (II. ip . (12. 3 p and (12. 4p . £/ien we /iawe /or any R> R 

EN A (X;R) = /3*EV A (X;R) (3.2) 
where E denotes the expectation with respect to Pp*. 

Proof. Using the GNZ formula ( 12. ip applied to the function h defined by (13.11) . we 
obtain 



EN A (X;R) = EJ2 Hpt\u)nB{u,R) = 0) 



ueX A 



= /9*E / t{Xn B(u,R) = 0)A(u,X)du. 

J A 

Now, using the finite range property (12. 3p . the assumption (12. 4p and the fact that 
R > R we continue with 

EN A (X;R) = /TE / 1(X n R) = 0)A(«, X n i2))du 

= P*E f l(Xn%fi)=|)%0)d« 
= i8*E7 A (X; 5). 

□ 



Using a classical ergodic theorem for spatial point processes obtained in 
we can expect that when the observation domain is large A^a(X; R) ~ E N A (X; R) 
and V A (X;R) ~ EV A (X.;R). Using this and Proposition 13. 1| we are suggested a 
natural and computationally cheap ratio-estimator of the parameter /3*. We define 
this estimator and present its asymptotic properties in the next section. 

3.2 Asymptotic properties 

To derive asymptotic properties of our procedure when the window expands to M. d , 
we assume that X is observed in A n where (A n ) n >i is a sequence of cubes growing 
up to R d . 

When u is a point close to the boundary of the window A n , we can only evaluate 
h(u, X n A n fl B(u, R)). One simple strategy for eliminating the edge effect bias is 
the border method, i.e. by considering the erosion of A n by R defined by 

A n,R = A n © R = {u € A„ : B(u, R) C A n }. 

Then we define our estimator by 

^ ^ „ N A „(X;R) 

P n = n (X; R) = "< R . (3.3) 

V K,r^ R > 

We let F denote the empty space function (see [14| ) given, in the stationary case, 

by 

F(r) = P^(XnS(O,r)^0) 
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for r > 0. We also define for u, v £ M. d the function F UjV : IR + — > [0, 1] by 

F u , v (r) = Pp* (X n B(u, r) ^ 0, X n r) ^ 0) 
for r > 0. We now state our main result. 

Proposition 3.2. Assume (ll.ip . (12. 3p and (12. 4p . and let R> R. Then we have the 
following statements as n — > oo. 

(i) The estimate n is a strongly consistent estimate of j3* . 

(ii) If Pp* is ergodic, the following convergence in distribution holds 

\K.~ R \ i/2 0n-n^m^ 2 ) (3.4) 

where 

a 2 = a 2 (r, R) = ^ ~ + f (1 - F 0;V (R))dv. (3.5) 



1 - F(R) (1 - F(R)) 2 J B (o,R) 
(Hi) The following convergence in distribution holds 

,1/2 0n ~n d 



\A~ R \ 1/2 ^—^^AT(0,l) (3-6) 



a 



n 



where a 2 , = ^(X; R) is a strongly consistent estimate of a 2 given by 

a 2 n = \A n] ,\ ^-^ + - KrK I ^ (3.7) 

' 1 W,JX:R) Vl_(X;R) ! 



with 

W A -{X;R) = [ [ l(Xn5(«,fl) =0)l(XnS(v,fl) = 0)dwdw. 

7A n jS JB(u,R)nA nJi 

If the measure P^* is not ergodic, then it can be represented as a mixture of 
ergodic measures (see [T0| Theorem 14.10]). In this case, the asymptotic distribution 
in (13. 4p becomes a mixture of Gaussian distributions. Since a mixture of standard 
Gaussian distribution is a standard Gaussian distribution, this explains why Pp* is 
assumed to be ergodic for the result (ii) and omitted for the result (iii) in Proposition 

GL2J 

We underline the interesting and unexpected fact that the asymptotic variance 
<r 2 (/3*,i?) (and its estimate) are independent on A which is the unspecified part of 
the Papangelou conditional intensity. 



4 Simulation study 

In this section, we investigate the efficiency of the Poisson intensity estimate in a 
simulation study. We consider the models (i)-(v) described in Section 12.21 We set 
R = 0.05 and 0* = 200 for all the examples except for the model g2 where we set 
0* = 50. We consider 
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• Strauss point process where models si and s2 have respectively 7 = 0.2 and 
7 = 0.8. 

• Strauss point process with hard-core with 5 = R/2 and where models shcl 
and she 2 have respectively 7 = 0.2 and 7 = 0.8. 

• Piecewise Strauss point process with p = 3 and (Ri, R2, R3) = (-R/3, 2R/3, R) 
and where models psl and ps2 have respectively (71,72,73) = (0.8,0.5,0.2) 
and (71,72,73) = (0.2,0.8,0.2). 

• Triplets point process where models tl and t2 have respectively 7 = 0.2 and 
7 = 0.8. 

• Geyer saturation point process with saturation threshold s = 1 where models 
gl and g2 have respectively 7 = 0.5 and 7 = 1.5. 

The estimate of the Poisson intensity we propose depends on a parameter R 
and ought to work for any value larger than R. To illustrate this, (3* is estimated 
for different values of R. We consider R = R x p, where p is a parameter varying 
from 0.9 to 1.2. Each replication has been generated on the domain [0,L] 2 and the 
estimates are computed on [R, L — R] 2 . To illustrate the convergence as the window 
grows up, the simulations are done with L = 1 and L = 2. The empirical results 
based on 500 replications are presented in Table [TJ 

The results are very satisfactory. We observe that for values of R > R (i.e. 
for values equal or larger than the finite range), the estimates are almost unbiased 
whereas a non negligible bias is observed when R < R. It is interesting to note that 
this bias is negative when the point process is repulsive and positive when the point 
process is attractive (i.e. the model g2). The closer R is to the finite range R, the 
lower the standard deviations are. The number of points with no i?-close neighbours 
is lower and lower as R grows up, which explains this variance increase. Finally, we 
can also highlight that the asymptotic results are partly checked. When the window 
varies from [0, l] 2 to [0, 2] 2 , the standard deviations of the estimates are divided by 
2 which turns out to be the square root of the ratios of the domains areas. 
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Mean (std dev.) for R = R x p 



Model 




Ti 
F 


= 0.9 




Ti — 1 


f 


= 1.1 


Ti 
F 


= 1.2 


si L 




1 


99 


174.0 


(26.2) 


203.6 


(34.2) 


204.3 


(38.3) 


205.5 


(44.4) 


L 




2 


393 


172.6 


(12.7) 


200.5 


(16.4) 


200.3 


18.5) 


200.6 


(20.8) 


s2 L 




1 


156 


192.0 


'31.3) 


203.3 


(41.8) 


203.9 


'48.1) 


204.7 


'55.8) 


L 


= 


2 


622 


190.7 


(16.0) 


200.7 


(19.5) 


200.3 


(22.8) 


200.4 


(26.5) 


shcl L 




1 


94 


174.3 


(23.2) 


201.9 


(30.3) 


203.3 


(36.8) 


205.0 


(41.9) 


L 




2 


379 


174.1 


(12.1) 


201.6 


(15.7) 


201.9 


(17.8) 


202.1 


'21.1) 


shc2 L 




1 


130 


194.3 


(30.2) 


205.4 


(36.7) 


207.0 


(42.8) 


209.1 


(51.6) 


L 


= 


2 


514 


190.2 


15.2) 


198.5 


(17.8) 


199.5 


(21.1) 


200.6 


(24.4) 


psl L 




1 


111 


172.2 


(25.1) 


202.6 


(33.2) 


203.5 


(38.7) 


206.5 


(47.2) 


L 




2 


445 


171.5 


(12.4) 


201.6 


(16.3) 


201.7 


(18.6) 


202.4 


(21.6) 


ps2 L 




1 


134 


191.1 


(30.7) 


201.7 


(37.7) 


206.2 


(45.2) 


208.8 


(54.2) 


L 





2 


535 


192.8 


(14.5) 


201.7 


(17.5) 


202.1 


(20.6) 


202.0 


(24.4) 


tl L 




1 


159 


181.6 


(34.6) 


204.5 


(40.6) 


205.2 


(48.7) 


205.8 


(58.8) 


L 




2 


634 


180.5 


(16.7) 


201.6 


(19.5) 


202.1 


(22.9) 


202.8 


(26.7) 


t2 L 




1 


178 


188.7 


(30.7) 


200.7 


(35.0) 


202.4 


(39.9) 


205.1 


(50.2) 


L 




2 


727 


189.3 


(17.5) 


200.0 


(20.2) 


199.8 


(24.5) 


199.8 


(30.8) 


gl L 




1 


110 


174.8 


(26.6) 


201.4 


(34.4) 


201.5 


(40.2) 


203.6 


(46.0) 


L 




2 


440 


175.0 


(13.3) 


200.7 


(17.4) 


201.3 


(19.9) 


201.6 


(23.1) 


g2 L 




1 


70 


54.5 


(10.4) 


50.0 


(10.4) 


50.2 


(11.4) 


50.5 


(12.1) 


L 




2 


276 


54.4 (5.4) 


50.2 (5.4) 


50.2 (5.8) 


50.3 (6.3) 



Table 1: Mean and standard deviations of estimates of the Poisson intensity pa- 
rameter /3* for different parameters R based on 500 replications of different mod- 
els generated on the window [0,L] 2 and estimated on the window [R, L — R} 2 for 
L = 1,2. The second column, n, indicates the average number of points (over the 500 
replications) falling into [0, L] 2 , i.e. the Monte-Carlo estimation of En(X fl [0, L] 2 ). 



A Proof of Proposition 13.2 



We use the following additional notation in the proofs. The notation stands for 
summation over distinct pairs of points, A c denotes the complementary set of A in 
M. d . Let r = (r x ) xeR d be the shift group, where t x is the translation by the vector 
— x G M. d . For brevity, when there is no ambiguity, we skip the dependence on X 



and on R of the variables N\ ~,V\ ~, W\ ~/ot. 

n,R n,R n,R ' l 



Proof, (i) There exists at least one stationary Gibbs measure. If this measure is 
unique, it is ergodic. Otherwise, it can be represented as a mixture of ergodic mea- 
sures (see [TU| Theorem 14.10]). Therefore, we can assume, for this proof, that Pp* 
is ergodic. We apply a general ergodic theorem for spatial point processes obtained 
by [15j. For this, wc have to check that N\ and V A arc additive, invariant by trans- 
lation and satisfy EAT Ao < oo and E V\ < oo, where for instance A = [0, l] d . The 
additivity property is straightforward. Regarding the invariance by translation, let 
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x G Q, A C R d and y G R d , we have 

N TyA (T y x;R) = t TyA (u)t(d(u,T y x\u) > R) 

= ~&a(u + y)t(d(u,T y x\u) > R) 

= J2 U(v)l(d(v - y, T y -K \ (v - y)) > R) 

= Mv)Hd(v, x\v)>R) = iV A (x; R). 

Following similar arguments, we also have V T \(t v x.] R) = V A (x; R). Let us now focus 
on the integrability conditions. First, we have 

E | Va I =EK Ao =E / l(Xn B(u,R) =0)d«=l- F(R) < oo. 

JAo 

Second, from Proposition 13. II 

E|iV Ao | =EN Ao = p*EV Ao = P*(l-F(R)). 
Therefore as n — > oo, we derive the following almost sure convergences 

Since the function (y,z) (-> | is continuous for any (y,z) in the open set {(y,z) G 

M 2 : z 7^ 0}, we deduce the expected result for any fixed R> R. 
(ii) From the definition of the innovations process (12.21) . we have 

, N A ~-p*V A ~ |A„s|- 1/2 / A -(X,/i) 
\K,n\ 1/2 ^ - = \K,n\ 1/2 \ " = u i-t^ • (A.2) 

"A g A nR 5 

Henceforth from Slutsky's Theorem and (lA.ip . the result will be deduced if we prove 
that |A n ^| -1 / 2 iA ~(X, ft,) tends to a zero mean Gaussian distribution. This will be 
done by combining two papers jl] and j3] dealing with asymptotic normality and 
variance calculation of innovations processes. For brevity, we write I A in place of 
/a(X». 

Following [4; Lemma 3], we have to assume that Pp* is ergodic and to check that 
for any bounded domain A: (I) E|/a| 3 < oo, (II) EIp — > for any sequence of 
bounded domains T n such that r n — > as n — » oo and (III) I A depends on X k ^. 
(Ill) follows from the finite range property (12.31) and the definition of h given by 
(13. ip . We now focus on verifying (I). We have from Cauchy-Schwarz inequality 

E|I A | 3 <4(EJVjJ + EV£). (A.3) 
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On the one hand, we have since for any m > 1, h m (u, x) = h(u, x) 

? 

EiV| =EiV A + 3E h(u,X\u)h(v,X\v) 

11,116X4 

+ E ^ h(u,X\u)h(v,X\v)h(w,X\w). (A.4) 

Let Ti , T2 and T3 denote the three terms of the right-hand side of (1A.4I) . Following 
the proof of (i) in Proposition 13. 2} ET X < /3*|A|. By using iterated versions of the 
GNZ formula ( 12. ip . we derive the following computations 

T 2 = 3E / / h(u,XUv)h(v,XUu)\p*({u,v},X)dudv 
JaJa 

= 3/T 2 E / / t(d(u, XUv)> R)t(d(v, X U u) > R)X(u, X)X(v, X U u)dudv 

J A J A 

= 3/T 2 E f [ l(d{u, XUv)> R)t(d(v, XUu)> R)\ 2 (u, 0)dwdt» 

J A J A 

< 3(3* 2 \A\ 2 . 

Using similar ideas, we leave the reader to check that T 3 < /3* 3 |A| 3 . On the other 
hand, we also check that 

EVI = [3* 3 e[ [ [ t(d{u,X) > R)t(d(v,X) > R)t(d(w,X) > R)dudvdw 

J A J A J A 

< /T 3 |A| 3 . 
We finally obtain 

E |/ A | 3 < 4/3*|A| + 12/3* 2 |A| 2 + 8/3* 3 |A| 3 < 00 

whereby (I) is checked. Following the proof of (I), we may prove that for any bounded 
domain A 

Ell < 2/3*|A| +4/3* 2 |A| 2 

which proves (II). Therefore we can invoke jU Lemma 3] to obtain the asymptotic 
normality of {3 n . 

The derivation of the asymptotic variance of |A n Var I\ g (X, h) is done using 
Lemma IB. II in Appendix IB1 with the function g = h given by (13. ip . We omit the 
details and leave the reader to check that h fulfills the integrability assumptions 
in Lemma fB.ll and focus on the derivations of the quantities A\,Ai and A3. First, 
it is clear that A\ = — F(R)). Second, note that from (12. 3p we have for any 
v E B c (0, R) and x E Q 

A^(0, x)\is*{v, x) - A^({0, v}, x) = 0. 

Since B(0,R) C B(0,R) 

E [ (\p,(0,X)\f ) *(v,X)-\p.({0,v},X))dv = 

J B(0,R)\B(0,R) 

11 



which means that we can substitute R by R in the term A 2 . Then, we derive 
A 2 =E [ l(d(0,X) > R)t(d(v,X) > R)(f3* 2 - X p *({0,v},X))dv. 

JB(0,R) 



Third, 



A 3 = E I (l(d(0, X U v) > R) - t(d(0, X) > R)) 

JB{0,R) 

x (l(d(u,XUO) > R)-t{d{v,X) > R))Xp*{{Q,v},X)dv 
= E [ l(d(0,X) > R)t{d{v,X) > R)\p*{{0,v},X)dv 



which leads to 



A 2 + A 3 = [3* 2 E f l(d(0, X) > R)t(d{v, X) > R)dv 

JB(Q,R) 

= /T 2 [ (l-F 0>v (R))dv. 
Jbio.r) 



In other words 



|A ni j|- 1 Var/ Ani5 (X,/ i )-^r(l-F(i?)) + r 2 / (1 - F 0jV (R))dv 

Jb(o,r) 



which leads to the result from (1A.1I) and (1A.2jl . 



(iii) We check that the proposed estimate o\ tends almost surely towards a 2 as 
n — > 00 (a convergence in probability would be sufficient). From (13. 7J) , (lA.lj) and 
(13. 5p . we need to prove that 



[ (l-F 0tV (R))dv. 
Jb(o,r) 



We decompose Wa - as W A _ + wl 2 ^ _ where 

n >R R n,R 



Wj? ] _ = [ [ t(d(u,X) > R)t(d{v,X) > R)dudv 
n ' H JA nJi eR Jb(u,r) 

/ / l(d(u, X) > R)t(d(v, X) > R)dudv 



where A n ^ i? = {u G A n ^ : B(u,R) C A n ^}. Regarding the first term, we can 
now apply the ergodic theorem to prove that almost surely 

\A nn eR\- 1 W { A 1] — ^E f l(d(0,X) > R)l(d(v,X) > R)dv 
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where the latter comes from the finite range assumption (12. 3p and the fact that 
R > R. Since |A _g| ~ \A n R i?| as n — >■ oo, the proof is done by verifying that 
I A k| _1 W / a 2 ' ) _ — > which is straightforward since 



\A ni ~ R \-^_ R < -* i |B( U>J R)| =o(l). 



□ 



B Auxiliary lemma 

The following result is a particular case of |5, Proposition 3.1]. 

Lemma B.l. Assume the Papangelou conditional intensity of the point process X 
satisfies ( II. ip and ( 12. 3 p finite range R. Let g : K d x — >■ R be a nonnegative 
measurable function satisfying for any u G M. d and xefi 

g(u, x) = 0(0, r„x) and g(u, x) = #(w, x B(u 5) ) 
with R> R. We also assume that the following integrability conditions are fulfilled 

E[0 2 (O,X)A j9 *(O,X)] <oo 

E / g(0, X)g(v, X) | A^(0, X)A^(^, X) - A^({0, v}, X)\dv < oo 
Jb(o,r) 



E / | (<?(0, X U «) - 0(0, X)) (<?(<;, X U 0) - </(«, X)) 
Jb(o,r) 



\p*({0,v},X)dv < oo. 



Then, as n — )■ oo 

IAJ" 1 Var J An (X, — > ^ + A 2 + A 3 

w/iere 

A 1 = E[/(0,X)A^(0,X)] 



A 2 = E / ^(0,X)(A^(0,X)A j a*(v,X)-A i a*({0,t;},X))d« 
Jb(o,r) 

A 3 = E [ (g(0, XUv)- g(0, X)) (<?(<;, X U 0) - g(v, X))A j9 *({0, v}, X)du. 
Jb<o,r) 



'B(0,R) 
1(0,R) 
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